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Abstract 

We investigate spherically symmetric equilibrium states of the Vlasov-Poisson system, 
relevant in galactic dynamics. We recast the equations into a regular three-dimensional 
system of autonomous first order ordinary differential equations on a region with compact 
closure. Based on a dynamical systems analysis we derive theorems that guarantee that 
the steady state solutions have finite mass and compact support. 
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1 Introduction 



In this paper we investigate spherically symmetric equilibrium states of the Vlasov-Poisson 
system. These equations describe a collisionless gas of particles that interact by the Newto- 
nian gravitational field they generate collectively. Examples of systems that can be modelled 
by these equations are ensembles of stars in e.g., a galaxy, a globular cluster or an ensemble 
of galaxies in a rich cluster 

A collisionless gas is characterized by the phase space mass density distribution f(t,x,v), 
where x, v £ M 3 and where / d 3 x d? v is the mass contained in the phase space volume element 
(i 3 x<i 3 v at (x, v). The spatial mass density is given by 

p(x) = J /(x,v)d 3 v. (1) 

The Newtonian potential V satisfies the Poisson equation AV = 4-7rp. Under the assumption 
that V tends to zero at infinity and with suitable fall-off and regularity conditions on the 
function p there is a unique potential V < for a given density p. It is assumed in the following 
that the function / is non-negative and spherically symmetric, i.e., that it is invariant under 
simultaneous rotations of x and v. The symmetry of / is inherited by p and V. Under these 
conditions the Poisson equation can be written as 

l(r 2 U'y = 4vrp(r), (2) 

where r = |x|, ' = d/dr, and where the gravitational constant has been set to 1. As the 
distribution function of a collisionless gas the function / satisfies the Vlasov equation. The 
following is concerned with time-independent solutions and in that case the Vlasov equation 
takes the form 

v-^ + VU-^ = 0. (3) 
ox ov 

Equations (|T|). (|2"|) and (j3J) constitute the Vlasov-Poisson system in the time- independent, 
spherically symmetric case. A definition which is useful when comparing models with colli- 
sionless matter and models in which the matter is described by a fluid is the radial pressure 
defined by 

J v 2 r f(x,v)d 3 v, (4) 
where v r = |v • x|/r is the radial velocity of a particle. 

For a time-independent spherically symmetric solution of the Vlasov-Poisson system the 
distribution function depends only on the energy per unit mass and particle E = ^|v| 2 + V(r) 
and the squared magnitude of the angular momentum L 2 = x x v| 2 . This statement was first 
proved in pQ under appropriate regularity assumptions 1 ; it belongs to a class of statements 
known as Jeans theorems in the astrophysical literature, cf. 0. Conversely, any function of 
E and L 2 satisfies the Vlasov equation. 

1 In PP and |1U| the function / was interpreted as a number density, but this leads to the same mathematical 
problem as studied here. 
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Now let R e (0, oo] denote the radius of support of the system, i.e., 



R = {inf r\ : p(r) = for r > n} (5) 

Because of physical applications, we are interested in systems in equilibrium that have finite 
total mass. It follows from Theorem 2.1 of jlUj that the monotonically increasing potential 
satisfies lim r ^R V(r) = Vr < oo and that there is a cut-off energy Eq = Vr such that 
f(E,L 2 ) = when E > Eq. We now introduce the definitions 

£ = E - E , uj = E - V = V R - V , (6) 

where u can be interpreted as a relative potential, and £ as the binding energy per unit mass 
for each particle. When f(£,L 2 ) ^ 0, then £ > 0, uj > 0, while £ and uj are zero at the 
boundary R of the system. 

In terms of £ and L 2 , the mass density takes the form: 



T 2 

"max 



p(r,uj) = ^l I f(£,L 2 )\v r \- L dL 2 d£ , (7) 



2tt 

'o Jo 

where L^ ax = 2r 2 (E - V(r)) = 2r 2 (w - £). 

The mass of the system inside the radius r is determined by 

m' = Anxr 2 p , (8) 

and thus m = 4ir J Q s 2 pds, assuming that there is no central point mass. As a consequence 
Poisson's equation takes the form {r 2 u)')' = m', which leads to 

-' = -^. (») 

Hence m and V are monotonically increasing functions of r while uj is monotonically decreas- 
ing. 

Given a suitable choice of the dependence of / on £ and L 2 the existence of corresponding 
global solutions with regular potential can be shown 1 . What is much harder is to decide 
whether these solutions satisfy the condition that the total mass is finite. One approach to 
this, which will be pursued in the following, is to try to prove the stronger statement that 
the radius of the support is finite. In ^H] finiteness of the radius was proved for a large class 
of solutions, which will be described in Section |21 Many models used in the astrophysical 
literature are covered by the results of jlUj . but some are not. Motivated in part by this we 
consider a wider class of models in the following and greatly extend the domain of validity 
of the results of ^Uj ■ The main assumptions used in ^H] concerned the behaviour of / near 
£ = 0. It turns out that if more general classes of solutions are to be handled it is often 
necessary to make assumptions concerning high values of £ as well. 

The mathematical problem to be solved is to obtain information on the qualitative behaviour 
of solutions of a system of ordinary differential equations. The technique which allows us to 
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go beyond what was done previously is to reformulate the problem using new variables, since 
this enables us to apply the theory of dynamical systems to the resulting equations. The 
geometrical intuition resulting from the dynamical systems formulation played an important 
role in the development of the proofs. The new variables are obtained by analogy with an 
approach of Heinzle and Uggla 5 to the Euler-Poisson system describing a self-gravitating 
fluid. 

The paper is organized as follows: in Section |2] we define and describe the class of distri- 
bution functions we study in this article. For these distribution functions, in Section |31 we 
reformulate the static Vlasov-Poisson system as a three-dimensional system of autonomous 
differential equations on a state space with compact closure. The system is subsequently 
analysed in Section |I] by using methods from the theory of dynamical systems; in particular, 
functions that are monotone along solutions of the dynamical system play a key role. Based 
on the results of Section 21 the main theorems are stated and proved in Section |SJ we formu- 
late conditions that guarantee finiteness of the radius R of solutions. We conclude the paper 
in Section El with some examples and remarks, and give an outlook on further applications of 
the techniques developed here. 



2 Distribution functions 

There exists a type of distribution functions that can be said to be among the mathematically 
simplest, the "generalized polytropes" 

/<*,*,={ - <f>°> , (1 o, 

where n, I and 0_ are constants with n > 1/2 and I > — 1. Under these conditions the 
mass density of these models is well-defined, p = p^r 21 uj n+l , where p_ = 2' +3//2 7r 3//2 r(^ + 
l)r(n - l/2)r(n + I + I)- 1 . The models (UUJ) exhibit invariance under scale transformations 
(of space and time), which corresponds to a symmetry of the Vlasov-Poisson equation; this 
leads to mathematical simplification. In particular, scale invariance gives rise to so-called 
homology invariants (e.g., [I], [S], [3]), i.e., quantities that are invariant under scalings, which 
are adapted to the symmetries of the problem. The use of homology invariants turns out to 
be crucial for our considerations. 

The mathematical analysis of static spherically symmetric solutions in the generalized poly- 
tropic case reduces to that of an ordinary differential equation known as the Emden-Fowler 
equation. This equation has a long history intertwining mathematics and astrophysics which 
we have not attempted to reconstruct. We note only that Fowler |Sj introduced this equation 
as a mathematical generalization of the Emden equation, the case I = 0, whose importance 
in the theory of stellar structure was well-known. Further early references can be found in [2] 
and pQ. It appears that the Emden-Fowler equation was introduced independently much 
later by Henon jo' : . , who derived it from a distribution function of the form introduced above 
and coined the term 'generalized polytropes'. 
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In this paper we discuss distribution functions that naturally generalize the polytropes: 

sv.V -{«■*>» %H , (H) 

where (/>(£) is a non-negative function which is measurable, bounded on compact subsets 
of the interval (0, oo), and integrable on [0,1]. Then, provided that I > —1, distribution 
functions of the type 1)11)1 give rise to a mass density 

p{r,u) = C l r 2l g l+1/2 {u), (12) 

where g m (m > —1) is defined as 

g m (u) := ^ <j>{£){u> - £) m d£ , (13) 
J o 

and where Q := 2 i+3 / 2 vr 3 / 2 r(/ + l)/r(/ + 3/2). The radial pressure of the models (p} is 
given by p ra d(r,uj) = Cir 2l g l+3/2 (uj)/(l + 3/2). 

It is natural to define a polytropic index function n(oS) according to 

n(oj ) = -l + — '- . 14 

a log a; 

For the special case of the generalized polytrope <j>{£) oc £ n ~ 3 / 2 the function gi+i/ 2 (u) is 
given by gi+i/zfo) oc uj n+l , and therefore n(u>) becomes a constant, n(uj) = n. 

We call distributions that lead to a (bounded) function n(uj) and satisfy n(uj) — > no for 
u — > asymptotically polytropic in the low ca; regime. Distribution functions of this type 
with no < 3 + I were treated in ^0j. One of the main aims of the present work is to treat 
more general distribution functions since this is needed for applications. 

If it is assumed in addition that 4> satisfies 4>(£) < const £ k for some k > —1 on a neigh- 
bourhood of £ = 0, then 5/ + i/2(^) is continuous for uj > 0. Moreover, when I > —1/2, 
9i+i/2( L0 ) S C 1 (0,oo), and thus n(u) is continuous for all uj > 0; compare with the results 
of ^Uj. For later purposes we require n(uS) to be of class C . In order to achieve this we 
assume that <j) is C 1 (0,oo) such that 4>'(£) is bounded on compact subsets of (0, oo) and 
< f > '{£) — const £ k for some k! > —2 on a neighbourhood of £ = 0. Then it can be shown that 
9i+i/2( UJ ) is C 2 for w > and hence that n(u) is C , provided that I > —1/2. For the proofs, 
and a discussion of the case — 1 < I < —1/2, see AppendixEJ The regularity conditions on 
6 stated above will be assumed from now on. 



3 Dynamical systems formulation 

The main idea of this work is to reformulate the static Vlasov-Poisson system, 

= 4nr 2 p(r, uj) (15a) 

dr 

-2 , irU 
— = — r m , (15b) 
dr 
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cf. 1(B)). Q, as a three-dimensional system of autonomous first order ordinary differential 
equations on a region with compact closure. 

The system (|15[l becomes an autonomous system of equations by regarding r as a supplemen- 
tary dependent variable and introducing a new independent variable As the next step 
we perform a transformation of variables from (m > 0, r > 0) to two dimensionless variables 

4Trr 3 p{r,uj) m 

u= , q= — , . 16 

m vuj 

Note that uq = 4irCir 2+21 gi + i/ 2 (o->) /o-> ■ Using the fact that I > — 1 it is possible to compute 
r in terms of u and q. Then m can be computed in terms of those variables and so the 
transformation is invertible for u > 0. 

Introducing £ := logr as a new independent variable and converting the system (|15[) to the 
new variables we obtain 

dill 

— = u (3 — u + 21 — n(u>)q — Iqj , (17a) 

g = q(-l + u + q) , (17b) 

duJ riv \ 

^ = , (17c) 

We proceed by defining bounded variables in order to obtain a dynamical system on a bounded 
state space. Recall that u, q, uj > and define bounded variables U, Q, (U, Q, Q) £ (0, l) 3 , 

1 + u 1 + q 1 + uJ V ' 

If a solution u(£), q((,), w(£) is given, a new independent variable A can be introduced by the 
relation d\/d^ = (1 — ?7) _1 (1 — Q)~ 1 . This yields a solution of the system of equations 

^ = C7(l - 17) [(1 - Q) (3 + 2/ - (4 + 2l)U) - (n(Q) + I) Q{\ - U)] (19a) 
cla 

^ = Q(l - Q) [(2Z7 - 1)(1 - Q) + Q(l - U)] (19b) 
^ = -0(1 - J2)Q(1 - [/) , (19c) 
where n(fi) is n(u)\^ n y 

The r.h.s. of the system ()19|) is C 1 in f2, because of the previous assumptions on the distri- 
bution function, and polynomial in U and Q. Hence, it is natural to smoothly extend the 
system to the side faces of the cube, so that the system of equations (|19jl forms a C 1 dynam- 
ical system on the state space [0, l] 2 x (0, 1). We observe that the side faces are invariant 
subspaces of the system, moreover, they contain the attracting sets for any orbit in the state 
space, as discussed in the next section. Note that a solution of (|19|) which is global in A need 
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neither correspond to a solution of (|17j) which is global in £, nor to a solution of (|15j) which 
is global in r > 0. 

For generalized polytropes n(f2) = const, and the equation for 0, decouples from the equations 
for U and Q. This is a direct consequence of the fact that the variables U, Q are homology 
invariants. However, in the case of a general distribution function, n{Q) gives rise to the 
three-dimensional coupled system (|19|) . 

4 Dynamical systems analysis 

Consider the dynamical system (|T9|) on the state space [0, 1] 2 x (0, 1). The side faces of the 
cube are invariant subspaces, where the fixed points of the system are located, cf. Figure ^ 
Tabled lists the fixed points together with the eigenvalues of the linearizations of the system 
at the fixed points. It is elementary to prove the facts summarized in this table. We observe 
that the flows on the side faces U = 0, U = 1, and Q = possess a simple structure. In each 
of these cases the evolution for the variable other than f2 does not contain £1 and the induced 
system does not depend on n(Q). 



n 




Figure 1: The state space [0, l] 2 x (0, 1) with fixed points. There exist four "fixed lines", i.e., 
lines consisting of fixed points (Li, L2, L3, L4); the position of L2 depends on /. 

Since the coefficients of the system are C 1 we can conclude that for any initial data there 
exists a unique local solution and that it extends to a global solution provided O remains 
bounded away from zero and one on any finite interval. The latter condition always holds 
since d/dX(logCl) and d/dX(log(l — fi)) are bounded along any solution. 

An important tool for the local analysis of dynamical systems near a stationary point is the 
reduction theorem (see jH], p. 48). It says that the flow of any C 1 dynamical system is 
topologically equivalent in a neighbourhood of a stationary point to a product of a standard 



7 



Fixed point 


u q n 


Eigenvalues 


Li 
L 2 
L 3 
U 


1 o n 

3+21 n q 
4+21 U U 

o n 

1 1 fl Q 


1,1,0 

3+2Z 1+2 n 
4+21 ' 2+1 ' U 

3 + 2/ , -1 , 

0,0,0 



Table 1: Local properties of the fixed points on the side faces. Here f^o is a parameter 
belonging to the interval (0, 1) in each case. 

saddle with the flow on any centre manifold. 

The lines L 2 and L3 are transversely hyperbolic saddles in the following sense. The lin- 
earization of the system about any point of one of these lines has a zero eigenvalue whose 
corresponding eigendirection points along the line. The other two eigenvalues are non-zero, 
real and have opposite signs. Each of these lines is a centre manifold for any of its points. By 
the reduction theorem this means that the flow near one of these points looks topologically 
like the product of a hyperbolic saddle with a line. The stable manifolds of points on L 2 and 
the stable and unstable manifolds of points on L3 lie in the boundary planes. It can be con- 
cluded from these facts that the boundary planes involved are mapped onto coordinate planes 
by the mapping which does the reduction. As a consequence any interior solution which has 
an w-limit point on L 2 also has w-limit points with Q = not lying on L 2 and w-limit points 
in the interior. Any interior solution which has an cj-limit point (a-limit point) on L3 has 
w-limit points (a-limit points) with Q = not lying on L3 and cj-limit points (a-limit points) 
with U = not lying on L3. We will show below in Proposition 14.11 that these properties of 
the transversely hyperbolic saddles in combination with the monotonicity properties of the 
system (j!9|l imply that an interior solution cannot have an w-limit point on L 2 , and neither 
an a- nor an w-limit point on L3. 

Similarly, L\ is a transversely hyperbolic source, where both non-zero eigenvalues of the 
linearization are positive and L\ is a centre manifold for any of its points. It follows from the 
reduction theorem that any point of L\ has a neighbourhood such that every interior solution 
which intersects that neighbourhood converges to a point of L\ as A — > —00 and no interior 
solution enters that neighbourhood from outside. 

Another important tool which we will use in the following is the monotonicity principle |llj . 
which will now be stated. Let <pt be the flow of a dynamical system on an open set U C K n 
with S C U an invariant set. Let Z be a C 1 function on S whose range is the interval (a, b) 
where a £ RU {— 00} and b € KU {+00} and a < b. If Z is decreasing on orbits in S, then for 
all x € S the a- and w-limit sets of the orbit starting at x are disjoint from S. Moreover it 
cannot be the case that Z tends to b at any w-limit point or that Z tends to a at any a-limit 
point. In 11^ this was stated only for U = W 1 but the same proof works for general U. 

Proposition 4.1. (Global dynamics). The uj -limit of every interior orbit is located on 
O = 0. 
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Proof. The proof of the theorem is based on the monotonicity principle. The function is a 
strictly monotonically decreasing function on (0, l) 3 and on the side faces U = (Q ^ 0) and 
Q = 1 [U 7^ 1). Accordingly, the monotonicity principle yields that the w-limit of an interior 
orbit must be located on U = 1, Q = or O = 0. Suppose that an interior solution has an 
w-limit point p with U = 1, Q > and S7 > 0. There is a solution on the boundary U = 1 
passing through p. Because of the simple structure of the dynamical system on the surface 
17 = 1 it follows that the new solution has an w-limit point with Q = 1 and 0, > 0. This 
implies that the original solution also has a limit point with XJ = 1, Q = 1 and O > 0. This 
limit point is on L4 and it will be shown in Lemma l4.2l that this is not possible. Suppose next 
that an interior solution has an w-limit point p with Q = 0, < U < 1 and £1 > 0. There 
is a solution on the boundary Q = passing through p. Because of the simple structure of 
the dynamical system on the surface Q = it follows that the new solution has an w-limit 
point on L2 with O, > 0. Hence the original solution also has an w-limit point with these 
properties. Because of the structure of L2 as a tranversely hyperbolic saddle this implies 
that the original solution has an interior w-limit point, in contradiction to what has already 
been proved. Suppose that an interior solution has an w-limit point on L3. Since L3 is a 
transversally hyperbolic saddle, the solution must also possess an w-limit point p with Q = 0, 
< U < 1; this case has already been excluded above. Finally, if an interior solution has an 
w-limit point on Li, then the fact that L\ is a transversely hyperbolic source also leads to a 
contradiction. Thus U = 1 and Q = cannot contain an w-limit of an interior orbit and this 
leaves f2 = as the only attracting set. ■ 

The fixed points on L4 constitute a special case because they are not only non- hyperbolic, but 
the linearization of the dynamical system ()19|) at such a fixed point has three zero eigenvalues. 

Lemma 4.2. No interior orbit has an a- or to-limit point on L4. 



Proof. We rewrite the system (|19p (where U < 1, Q < 1) in cylindrical coordinates centred 
at (1,1,0). Define 



The coordinate <fi ranges in (0,7r/2), f is in (0,f max ), where < f max < 1 can be chosen 



f = y/(l - U) 2 + (1 - Q) 2 4> = arctan 



1-Q 



z = n. 



(20) 



1 - U 




The fi(f, (ft, z) are of the form 



f r (f, (ft, Z) = f r ,l{(ft, Z)r + f r ,2{4>, z)r 2 + / r , 3 (^, z)f 3 

U(r, (ft, z) = Ufi{(ft, z) + z)f + /^ i2 (0, z)f 2 



(22a) 
(22b) 
(22c) 
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where the fij are polynomials in cos eft, sin eft. Thus the system can be smoothly extended 
to also include the boundaries f = 0, (j> = 0, and <fi = tt/2, so that the state space is 
[0, f max ) x [0, 7r/2] x (0,1). From the construction it follows that the subset f = can be 
regarded as a blow-up of the fixed line L4; clearly, f = is an invariant subset. 

A fixed point analysis of the system reveals that there exist only the fixed points (0, tt/2, zq) 
(zq € (0, 1)) on f = 0; elsewhere z is monotone. These fixed points turn out to be transversely 
hyperbolic saddle points and hence no orbit with r > can have an a- or w-limit point with 
r = 0. This establishes the claim of the lemma. ■ 

Proposition 4.3. Assume that n(Q) < 3 + 1 for all fi < Qq (for some arbitrary < Qq < 1). 
Then the u -limit of every interior orbit lies on Q = 1, $7 = 0. If —I + e < n(Q) < 3 + I for 
all £1 < f^o (for some e > 0) then the uj -limit of every interior orbit is the point (0,1,0); if 
n(Q) < —I — e for all fi < (for some e > 0), then the co-limit is (1, 1,0). 



Proof. Consider the function 



The function Z is strictly monotonically increasing on all interior orbits, 

— = (2(1 + l)U(l - Q) + (3 + 1 - n)Q(l - U))Z > . (24) 

Since dlogZ/dX > it follows that Z tends to a limit, finite or infinite as A — > 00. In the 
latter case, Z — > 00, we observe (1 — U)(l — Q) — ► 0, so that the w-limit of the orbit must lie 
on (U = 1) U (Q = 1). Now assume the first case, lim^-^oo log Z < 00. From (|2*1)) we obtain 
dlog Z/dX > 2(1 + l)U(l - Q), hence 

/oo 
17(1 - Q)dX < 00 . (25) 

Since, firstly, £7(1 — Q) > and, secondly, the derivative of U(l — Q) is always bounded, we 
conclude that U(l — Q) — > as A — > 00. Assume that Q /> 1. Then there is sequence X n 
such that 1 — Q(X n ) > const > and U(X n ) — > for n — > 00. This implies that Z(X n ) — > 0, 
a contradiction. Thus, Q — > 1 in the limit A — > 00, and the w-limit of the orbit must lie on 
Q = l. 

In combination with Proposition 14.11 this yields that the w-limit of every interior orbit must 
lie on (U = 1) n (O = 0) or (Q = 1) n (O = 0). 

Consider the set S 1 = {(U,Q,tt) \ (2U - 1)(1 - Q) + Q(l - U) > 0} in the state space, see 



Figure 2(a) Note that U = 1 and Q = 1 are contained in Si. It is not difficult to show that 
Si is a future invariant set: firstly, note that dQ/dX = on the boundary of Si, cf. (|19b|) . 
Secondly, we observe that 

= U(l-U)(l-Q)[3 + l-n + 2(n-2)U}>0 on dSi , (26) 

cLX 
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so that dS\ acts as a "semipermeable membrane" . 

On S*i the function Q is strictly monotonically increasing except when Q = and Q = 1. 
Application of the monotonicity principle yields that the u;-limit of every orbit must lie on 
Q = 1. Hence the first part of the proposition is established. 

Assume — / + e < < 3 + / on (0, fio] (for some e > 0). Consider the set S2 = 

{(U, Q,£l) \Q > max sup^ {z+2i)+\i+n) ) } ^ n t^' x ^°) - ^ 2 * s future invariant, which 
is because d£l/d\ < and dQ/dX > on 95*2. On 5*2 the function U is strictly monotonically 
decreasing except when U = or U = 1. Accordingly, the monotonicity principle yields that 
the w-limit of every orbit must lie on U = 0. Combining this result with the previous state- 
ment we see that the o;-limit of every interior orbit is located on (Q = 1) D = 0) D (U = 0), 
i.e., it is the point (0, 1,0). 

Assume n(S7) + / < — e on (0, ^o] (for some e > 0) and consider the set S3 depicted in 
Subfigure |2(b)| Its boundaries are Q = 1 — 5 (for some appropriate 5) and the surface given 
by 

U = (3 + 2Q(l-Q)-eQ 

(3 + 2/)(l-Q)-eQ+(l-Q) ' 1 ' 

S3 can be shown to be future invariant, and dU/d\ > on S3. Using the monotonicity 
principle in combination with the previous statement of the proposition, we see that the 
w-limit on every orbit must be the point (1, 1,0), as claimed. ■ 

Proposition 4.4. Assume that n(£l) < 5 + 31 for all < Oo (f or some < f^o < 1 ) and 

n{Q) ^5 + 3/ in a neighbourhood of O, = 0. Then the co-limit of every orbit originating from 
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Li n (Q < n ) or L 2 n (O < n ) lies on Q = 1, O = 0. If -I + e < n(Q) < 5 + 3/ f/or some 
e > /or a// small enough, then the co-limit of every orbit is the point (0, 1,0). 



Proof. Consider the function 

JJ \ 2(1^77 



1 - u 



1 



Q 



l-Q 



3+21 
2(1+0 



1 



Q 



u 



l-Q (3 + 2Z)(1 - Z7) 



(28) 



denned for (U, Q) G [0, 1) x [0, 1) and consider the surface <J?(/7, Q) = in the state space. In 
the case n(f2) = 5 + 3/ the function <I> is a conserved quantity of the flow. Under the given 
assumptions the set (<3? > 0) n (S7 < f^o) is a future invariant subspace of the state space 
because ($ = 0) H (O < f^o) is a semipermeable membrane: 













d$ 






n{Q) 


~~ dX 


n(n) 


" ~dX 




a$ 


fdU 




dU 






~ dU 


\~dX 




~ ~d~X 


5+31 J 



n=5+3l 



+ 





fdQ 




dQ 




dQ 


[dX 




dX 


5+3l) 



Now, 



<9$ 

dU 



m* 1 



Q 

4 



U) 2 Q [(5 + 3/) -n(fi)' 



2(1+0 



( 2^-l)(l 



Q) + Q(l-U) 

QY{i-uyu 



(29) 



(30) 



(l-EO(l-Q) 

From dSOJ) we see that d$/dU > if and only if (U, Q) € Si, cf. Fig. |2(a)| Since $ = is 
contained in Si, d^/dU > and therefore d<&/dX > 0, as claimed. The surface $ = and 
the set $ > are depicted in Figure |31 




Figure 3: The future invariant set ^ > for 1 = 1. 

Since the set (<£> > 0)n(Q < Qq) is future invariant, every orbit originating from LiC\(Q < Qq) 
in entirely contained in this set. The same is true for orbits that originate from L 2 but this 
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is more delicate to prove, since L2 lies on the boundary of the set > 0. Recall, first, that 
dQ/dU is positive on S\, which is a neighbourhood of the surface (<3? = 0)\{(0, 1/2, fi) | f2 € 
(0, S7q]}. Now consider a solution that starts from L2 with 0, = f2i < ^q- Along the solution, 
$ — > for A — > —00 holds. The orbit is in Si as long as it is sufficiently close to L2, i.e., for 
sufficiently large negative A we observe dQ/dU > 0. If Vt\ is such that n(VLi) < 5 + 3/, then 
d^/dX > in the limit A — > —00, hence d<&/dX > for all A in a neighbourhood of —00. We 
conclude that $ is positive for all large negative A and thus for all A, since > is a future 
invariant set. If n(Oi) = 5 + 3/ but n(fl) < 5 + 3/ for any £1 < Sli then the same argument 
applies. If neither of these conditions are satisfied let be the smallest number for which 
n(£l) = 5 + 31 on the interval [f^^i] and let A2 be such that ^(^2) = ^2- For A < A2 the 
solution lies on the surface $ = 0; however, d^/dX > for A slightly larger than A2, since 
the solution is in S%. Thus in all cases <& eventually becomes positive. 

To continue the proof we proceed as in Proposition 14.31 Since (<& > 0) C Si, Q is strictly 
monotonically increasing along interior solutions in (<]? > 0) H (fl < Oo); by the monotonicity 
principle the u;-limit of an interior solution must lie on Q = 1. This proves the first assertion 
of the proposition. 

To conclude the proof, recall that U is strictly monotonically decreasing on S2 when n(p,)+l > 
e holds, which is the case for all small f2 by assumption. The monotonicity principle yields 
that the w-limit of an interior orbit is the point (0, 1,0). ■ 

Proposition 4.5. ( Global dynamics) . The a-limit of an interior orbit is a fixed point on 
L\, L2, or contained in Q = 1. 

Proof. The function fl is a strictly monotonically decreasing function on (0, l) 3 , and on the 
side faces U = (Q / 0) and Q = 1 (U ^ 1). Accordingly, the monotonicity principle yields 
that the a-limit of an interior orbit must be located on U = 1, Q = 0, or Q = 1. The line 
L4 cannot act as an a-limit set for interior orbits, see Lemma 14.21 Suppose that an interior 
solution has an a-limit point p with U = l, 0<Q<1 and Q < 1. Then the entire orbit 
through p must be contained in the a-limit set, and hence a point on L4 because of the 
structure of the flow; a contradiction to what has already been shown. Therefore any a-limit 
point with £1 7^ 1 satisfies Q = 0. Suppose that an interior orbit has an a-limit point on L3. 
Since L3 is a transversally hyperbolic saddle, the orbit must also possess an a-limit point 
with U = 0, Q > 0, a contradiction. On Q = 0, a-limit points with (3 + 21) /{A + 21) ^ U < 1 
are excluded, since L\ is a tranversely hyperbolic source and because of the fact that the 
transversely hyperbolic saddle L3 cannot act as an a-limit set as shown above. Therefore, 
any a-limit points with 0^1 must be on L\ or L2. ■ 

If a solution of the Vlasov-Poisson system possesses a regular potential V, i.e., if V(r) is 
C 1 on [0,oo), then Q, converges to some limit fio G (0, 00) (and thus to — > ujq) as r — > 0, 
whereby Q -> 0, and U -> (3 + 21) /(A + 21) (r -> 0). To show this we use (Q1J and 
fi7+i/2( w ) = 52+1/2(^0) + o(l) as r — > 0. The corresponding solution (U, Q,0)(A) of (fl9|) 
exists globally to the past and the above limits are attained as A — ► —00; the solution is 
therefore associated with an orbit whose a-limit resides on L2. Although they have a regular 
potential, these solutions do not necessarily possess a regular density function p(r) since 
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p oc r 2l gi +1 / 2 (^o) as r — > 0, and therefore /? — > oo for Z < 0. Nevertheless, for reasons of 
brevity, we refer to these solutions as regular solutions. 



5 Finiteness of the radius of the support 

For the theorems, let us restate the regularity requirements on <f>: we assume that 4>(£) is 
C 1 (0, oo), and that cf), cj)' are functions, which are bounded on compact subsets of (0, oo). It 
is required that 4>(S) < const £ k for some k > — 1 and <p'(£) < const £ k for some k' > —2 on 
a neighbourhood of £ = 0. When I < —1/2, £<f>(£) must be Holder-continuous with an index 
strictly greater than — I — 1/2. 

Theorem 5.1. Let <f>(£) be such that n(u) < 3 + I for all uj < ujq (for an arbitrarily small 
(Jo > 0). Then every associated solution of the Vlasov-Poisson system has finite mass and 
radius. 

Proof. We investigate the relation dr/dX = r(l — U)(l — Q). Integration yields 



R = r exp 



(l-U(X))(l-Q(X))d\ 



U 



(31) 



where Ao, r$ are such that (U, Q, 0)(Aq) corresponds to (m(ro), w(ro), ro). Thus, to prove 
that R < oo, it suffices to show that the integral in (|31j) is finite. 

From Proposition 14.31 we know that the w-limit of every interior orbit is a subset of the set 
(O = 0) fl (Q = 1). Define 5Q = 1 — Q and 5U = 1 — U. Since Q is strictly monotonically 
increasing on the set Si, see Figure 2(a)[ 5Q is strictly monotonically decreasing on every 



orbit for sufficiently large A. Written in 5Q, SU, Eq. (|19b|) reads 

&1 = -5Q [SU(1 - 4£Q) + 6Q(1 - 5Q + 35U5Q) } . (32) 
dX 

Since 5Q — * for A — > oo, for every small e > there exists some A £ such that 

d{6Q ) 
dX 

for all A > X F . 



< -5Q [(1 - e)5U\ (33) 



Consider the integral 5U5Q dX appearing in (|31|). where Ao is chosen to be greater than 
A e , where e is small. We obtain 

oo <5Q(A ) _ x SQ(Xq) 

J 5U5QdX= J 5U5Q (-^^J d(5Q) < J (1 - eyHtfQ) < oo , (34) 

A 

i.e., the integral is finite. ■ 
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Theorem 5.2. Let (f>(£) be such that n(uj) < 5 + 3/ for all small uj, and define oo CT \t '■= 
sup{w | n{u) < 5 + 3/} ; and w c := Vr — V(0). Then every associated regular solution with 
w c < w cr it has finite mass and radius. 

Proof. Recall that for a regular solution the a-limit set must be a point on L<i- With this 
information, the statement of the theorem can be deduced from Proposition 14.41 in the same 
way that Theorem 15,11 was deduced from Proposition 14.31 ■ 

6 Conclusions and outlook 

In this paper we have investigated stationary solutions of the Vlasov-Poisson system. We have 
derived conditions on the distribution function guaranteeing that the resulting steady states 
have compact support and thus finite total masses. Theorem 15.11 formulates a condition 
that is local in the sense that only the properties of the distribution function for particle 
energies close to the cut-off energy are relevant; the theorem is a generalization of a theorem 
proved in ^U]. The condition of Theorem l5.1l is not only sufficient but also necessary 2 for the 
statement to hold. The condition stated in Theorem 15. 21 is not local; this criterion is based on 
the properties of the distribution function for a finite range of particle energies. Theorem 15.21 
thus opens access to a class of distribution functions that could not be treated previously. 
The condition of the theorem is sufficient, but not necessary. 

In the astrophysical literature there exists a large variety of stationary spherically symmetric 
models built of self-gravitating collisionless matter. Many of these models are based on 
distribution functions for which the theorems of this paper are relevant. As examples we 
mention the King models, the Wooley-Dickens models, see |2|, and the Kent-Gunn models, 
see [JJ: Theorem 15. II applies. 

Distribution functions of the form (fi(£) oc (exp£ — 1 — £)L 21 — the isotropic case is called 
the Wilson model [2] — are in general not covered by Theorem 15.11 However, Theorem 15.21 
applies for a wide range of /, since n(uj) — > 7/2 as to — > 0. Theorem 15.21 ensures finiteness 
of the radii and masses of solutions that satisfy lo c < u> C rit; since n(u>) is an increasing 
function which diverges asw-> oo, to CT \t is necessarily a finite number. The theorem does not 
give information about radii and masses of regular solutions with uj c > w C rit- A numerical 
investigation suggests the following behaviour: when i = all regular solutions, including 
those with u> c > uJcrit, possess compact support, which is in agreement with [HE} . However, for 
a certain range of values of I, the situation is more complex: for almost every value uo c > w cr it 
the associated regular solution has a finite radius, but there exists a discrete set of values 
{uj Ct i | i = 1 . . . n, 1 < n < oo} such that the associated solutions extend to infinity. Hence, in 
many cases, but not in general, the failure of the condition n{uj) < 5 + 31 to hold for large u> 
entails the occurence of regular solutions that are infinitely extended. 

2 When we restrict the theorem to the class of asymptotically polytropic distribution functions, this follows in 
a straightforward manner from a local dynamical systems analysis; we refer to the paper where analogous 
issues occured in the context of static perfect fluid solutions. Generalizations to other classes of 4>(£) are 
possible, but we refrain from discussing them here. 
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The above examples of distribution functions are all asymptotically polytropic in the low 
u regime; this allows one to also include the boundary $7 = in the dynamical systems 
state space and to obtain specific information about the asymptotic properties of solutions. 
When the models are asymptotically polytropic or "asymptotically isothermal" in the high 
ijj regime, it is likewise possible to include Q = 1 (although the isothermal case requires a 
slight change of variables). In this context, one might be able to derive additional theorems 
concerning mass-radius properties of solutions in analogy to what was done in the perfect 
fluid case in [3], e.g., it might be possible to establish mass-radius relationships. 

In the isotropic case the static Vlasov-Poisson system coincides with the static Euler-Poisson 
system which describes equilibrium states of self-gravitating perfect fluid matter. The colli- 
sionless matter model with density p and radial pressure p r ad> given in Q and ©, can thus 
be interpreted as a perfect fluid matter model with density p and pressure p = p ra a with 
a barotropic equation of state p(p) implicitly determined by the distribution function <j). It 
is straightforward to construct theorems that are analogous to 15.11 15.21 for the perfect fluid 
case; such theorems cover more general equations of state than the class of asymptotically 
equations of state discussed in [3]. 

In the present paper we have treated a collisionless gas in Newtonian theory of gravity. 
However, the methods used here are likely to be of relevance also in the general relativistic 
case; the relativistic theorems in jlUj might thus be generalized. 



A Assumptions and properties of (f)(£), g m (cu), and n(co) 

In this section we investigate the regularity assumptions on <j) that guarantee that n{uo) 
becomes a (^-function for to > 0. 

By setting x = 6 /to in the definition (J13|) of g m (uj), we obtain 



i 



g m (uj) = u m+l / 4>{ux){l - x) m dx (m > -1) . (35) 



o 



First some identities are obtained under the assumption that <j) is smooth. Then rougher 
functions cf> are treated by approximating them with smooth ones and passing to the limit. 
If cj) is smooth then differentiating 1)35|) implies 



_d_ 

duj 



g m (u) = (m + l)uj m J <P{ujx){1 - x) m dx + iv m+1 J x</>'(ux)(l ~ x) m dx. (36) 



o 

Integrating by parts gives the identity 
d 



^9m{^) = uj m <p{ujx)x(l - x) m + mg m ^i(uj) = mg m _ l (uj) , (37) 



o 



for m > 0. For m = the result dgo{oj)/duo = 4>{oo) is obtained. 
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The case —1 < m < is more subtle. Each of the integrals on the right hand side of (|36|) 
can be split into the sum of an integral from zero to 1 — e and an integral from 1 — e to one 
for any e between zero and one. As e — > the second term in each of the integrals tends to 
zero. The first term in each of the integrals can be treated as was done for the integral from 
zero to one above. As a result 



_d_ 



u m lim 



co m lim 



co m lim 



1-e 



1-6 



(ujx)x{1 - x) m ^ +m J (f)(uJx){l-x) m - L dx 



1-e 



e m (Muj-euj)-6(uj)) +e m (b(u 



) + m J 4>{ujx){ 




1 - x) m - l dx 



l-e 



[uj — euj) — 0[u) 



+ 4>{uj)-m J [4>(uj)-4)(ujx) ] j{l-x) m - 1 dx 


e5 <Hu - eu) - _ m y - fax) _ xrl+S(Lr 



co m (l)(io) + uj m lim 



(1 - x) a 



where a is defined as a = —m + 5 for an arbitrary 5 > (but preferably small). 

Now suppose that <f> is a measurable function, which is bounded on compact subsets of (0, oo) 
and satisfies <fi(£) < const £ k for some k > —1 on a neighbourhood of £ = 0. Then £~ k (j) is 
bounded on compact sets. It is possible to approximate it on any compact subset of [0, oo) by 
a sequence of smooth functions which are uniformly bounded and converge to it pointwise. 
By the dominated convergence theorem g m {u) is continuous and the sequence of functions 
g m defined by the approximants converge pointwise to that defined by <f>. Since (|H7|) holds 
for each of the approximants it also holds for cj) itself. It can be concluded that g m is C 1 . In 
the case m = 0, if is continuous then g m is C . Equation (|H6j) can alternatively be written 
as 



d_ 

dw^ 



l (uj)=LU- 1 ((m+l)g m (u)+LU m+1 I {ujx)(j)' {ujx){1 - x) m dx \ 



(38) 



If (j) is C 1 then we can use 1)38(1 for each of the approximants and pass to the limit to get the 
same relation for <p. Accordingly, the function g m becomes C 2 if the second term in 1)38(1 is 
continously differentiable. We note that this term is of the form 1)35(1 where (f>(y) has been 
replaced by y4>'{y). Using this formal resemblance we can simply adapt the conditions on (j) 
discussed above; hence, g m (u)) € C 2 (0, oo), if £<ft'(£) is bounded on compact subsets of (0, oo) 
and £(/)'(£) < const £ k for some k > —1 on a neighbourhood of £ = 0. 

Now consider the case m < 0. When we assume that there exists 5 > such that <p G C a , 
a = —m + 5, i.e., such that is Holder-continuous with index a, then the approximants for 
(ft can be chosen to converge in the appropriate Holder norm and we can use the identity for 
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the approximants and pass to the limit to get the corresponding relation for (p. Then 



and it follows that g m is C 1 . 
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